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Abstract 

In this paper we describe a measurement technique capable of resolving femtosec- 
ond X-ray pulses from XFEL facilities. Since these ultrashort pulses are themselves 
the shortest event available, our measurement strategy is to let the X-ray pulse sam- 
ple itself. Our method relies on the application of a "fresh" bunch technique, which 
allows for the production of a seeded X-ray pulse with a variable delay between 
seed and electron bunch. The shot-to-shot averaged energy per pulse is recorded. It 
turns out that one actually measures the autocorrelation function of the X-ray pulse, 
which is related in a simple way to the actual pulse width. For implementation of 
the proposed technique, it is sufficient to substitute a single undulator segment 
with a short magnetic chicane. The focusing system of the undulator remains un- 
touched, and the installation does not perturb the baseline mode of operation. We 
present a feasibility study and we make exemplifications with typical parameters 
of an X-ray FEL. 



1 Introduction and method 

The measurement of X-ray pulses on the femtosecond time scale constitutes 
an unresolved problem. It is possible to create sub-ten femtosecond X-ray 
pulses from XFELs, but not to measure them. In fact, conventional photode- 
tectors and streak-camera detectors do not have a fast enough response time 
to characterize ultrashort radiation pulses. For example, the rise time of the 
best streak-cameras approaches 100 fs, far too slow to resolve femtosecond 
pulses. Special measurement techniques are needed. In this paper we pro- 
pose a new method for the measurement of the duration of femtosecond 
X-ray pulses from XFELs. The method is based on the measurement of the 
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scheme for characterization of femtosecond X-ray pulse from the intensity autocorrelation measurements 
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Fig. 1. Experimental layout for ultra-short X-ray pulse-measurement using "fresh" 
bunch technique. Pulse and electron bunch are separated at the exit of first part 
of undulator, one is variably delayed with respect to the other. Modulation of 
the electron bunch is washed out. X-ray pulse and "fresh" electron bunch are 
overlapped in the second part of undulator. The averaged X-ray pulse energy at 
the exit of setup is measured versus delay yielding the autocorrelation trace. 

autocorrelation function of the X-ray pulses. The setup in Fig. 1, similar 
to that described in fl]], may be used to this purpose. The electron bunch 
enters the first part of the baseline undulator and produces SASE radiation 
with ten M W-level power. After the first part of the undulator, the electron 
bunch is guided through a short magnetic chicane whose function is both, 
to wash out the electron bunch modulation, and to create the necessary 
offset to install an X-ray optical delay line. The chicane is short enough to 
be installed in the space of a single XFEL segment, as shown in Fig. 2, and 
does not perturb the focusing structure of the machine. The optical delay 
line is sketched in Fig. 3, and has already been discussed in ||2|]. 

After the chicane, the electron beam and the X-ray radiation pulse produced 
in the first part of the undulator enter the second part of undulator, which 
is resonant at the same wavelength. In the second part of the undulator the 
first X-ray pulse acts as a seed and overlaps to the lasing part of electron 
bunch. Therefore, the output power rapidly grows up to the GW-level. First 
and second undulator parts are identical and operate in the linear FEL 
amplification regime. The relative delay between electron bunch and seed 
X-ray pulse can be varied by the X-ray optical delay line installed within the 
magnetic chicane, as illustrated in Fig. 4. Within the ID FEL theory, which 
is not too far from reality for the SASE X-ray case with a large diffraction 
parameter, one can write the shot-to-shot averaged power in the pulse from 
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Fig. 2. Installation of a magnetic delay in the baseline XFEL undulator. 
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Fig. 3. X-ray optical system for delaying the SASE pulse with respect to the electron 
bunch. The X-ray optical system can be installed within the magnetic chicane 
between first and second part of the undulator. 

the first part of the undulator as 



<P(0) = P exp[2LJRe(A(i))] (1) 

where P is the equivalent shot-nose power, L w = 36 m (6 cells) is the length 
of the two identical undulator parts and Re(A(f)) is the time-dependenl[£l 
field growth-rate. Similarly, the shot-to-shot averaged power in the pulse 
from the second part of the undulator, which is seeded with (P(t - t)), t 



2 Here we deal with a parametric amplifier, where the properties of the active 
medium, i.e. the electron beam, depend on time. 
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Fig. 4. Sketch of X-ray intensity autocorrelation measurement using tunable X-ray 
optical delay line within magnetic chicane. This autocorrelator produces complete 
autocorrelation trace. 

being the variable delay, can be written as 



<P 2 (f, t)> = (P(t - t)> exp [2L w Re(A(f))] = 1 <P(f - x)> <P(f)> . (2) 

"o 

The subsequent measurement procedure consists in recording the shot- 
to-shot averaged energy per pulse at the exit of the second part of the 
undulator as a function of the relative delay between electron bunch and 
seed X-ray pulse, with the help of an integrating photodetector. This yields 
the autocorrelation function 



A(t) = J dt (P(t — t)) (P(f)) . (3) 

— oo 

Autocorrelation measurements are well known methods in laser physics. 
Early on, it was realized that the only event fast enough to measure an ul- 
trashort pulse is the pulse itself. A number of schemes have been developed 
over the past decades to better measure ultrashort laser pulses. Most of them 
have been experimental implementations and variations of autocorrelators, 
i.e. devices capable of measuring the autocorrelation function of a given 
pulse, Eq. (3). Our scheme actually provides a device capable of performing 
an intensity autocorrelation measurement. 

Note that, in order to perform intensity autocorrelation measurements, one 
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Fig. 5. Experimental layout for measuring the optical pulse intensity versus time. At 
optical frequencies, an intensity autocorrelator using second-harmonic generation 
can be used. A pulse is split into two, one is variably delayed with respect to 
other, and the two pulses are overlapped in a second harmonic crystal. The second 
harmonic pulse energy is measured versus the delay, giving the autocorrelation 
trace. 

should insert a nonlinear element into an interferometer. In ultrashort laser 
physics, the most common approach in the visible range involves second- 
harmonic generation (SHG), in which a nonlinear crystal is used to generate 
light at twice the input optical frequency (see Fig. 5). The measurement 
procedure is to record the time-averaged second harmonic pulse energy 
as a function of the relative delay % between the two identical versions 
of the input pulse. Due to nonlinearity, the total energy in the second- 
harmonic pulse is greater when the two pulses incident on the nonlinear 
crystal overlap in time. Therefore, the peak in the second-harmonic power 
plotted as a function of % contains information about the pulse width. 

Similarly, here we assumed that we dispose of an ensemble of identical 
pulses, so that the autocorrelation function can be constructed from a large 
number of energy measurements taken for a different delay parameter t. 
The measured energy is the sum of a constant background term due to 
startup from shot noise from each part of the undulator, and of the intensity 
autocorrelation term, which arises from the interaction of the delayed elec- 
tron bunch and the seed SASE pulse from the first part of undulator. Due 
to high-gain FEL amplification in the second part of the undulator, the total 
energy in the X-ray pulse at the exit of the setup is much higher when the 
seed SASE pulse and the electron bunch overlap in time. This means that we 
effectively deal with a background-free intensity autocorrelation function 
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measurement. Therefore, the peak in the shot-to-shot averaged energy of the 
X-ray pulse at the setup exit, plotted as a function of t, contains information 
about the averaged X-ray pulse width. 

One immediately recognizes the physical meaning of the autocorrelation 
function. The Fourier transform of the autocorrelation function, A(co), is 
related to the Fourier transform of the signal function W(co), i.e. to the Fourier 
Transform of the intensity vs time, by A(co) = \W(co)\ 2 . An autocorrelation 
function is always a symmetric function. Thus, A(co) is a real function, 
consistent with a symmetric function in the time domain. The intensity 
autocorrelation function assumes its maximum value at t = 0. Moreover, 
the autocorrelation function is an even function of x, independently of the 
symmetry of the actual pulse. Therefore, one cannot uniquely recover the 
pulse intensity profile from the knowledge of the autocorrelation function 
only. This is also understandable from the fact that correlation techniques 
provide the possibility to measure the modulus of the Fourier transform of 
the signal function, while information about its phase is missing. 

However, since the pulses exhibit no overlap for delays much longer than 
the pulse width, the autocorrelation function goes to zero for values of t 
larger than the pulse width. Therefore, the width of the correlation peak 
gives information about the pulse width. One can estimate the FWHM of 
the radiation pulse from the knowledge of the autocorrelation function if 
one assumes a specific pulse shape. Then, the FWHM can be found by divid- 
ing the intensity autocorrelation FWHM by a deconvolution factor, which is 
specific for a given shape. If one deals with smoothly varying pulse shapes, 
the deconvolution factor is about 1.5 |0r|8t] and the variation in the decon- 
volution factor is of the order of 10% only. Therefore, the pulse duration 
can be approximatively obtained from the knowledge of the FWHM of the 
intensity autocorrelation function, even though the pulse shape remains 
unknown. 

In the next Section we will present a feasibility study of our method, based 
on simulations with the code Genesis 1.3 [9], with the help of parameters in 
Table 1. We further discuss an alternative method for radiation pulse width 
measurement which is based on simpler hardware, but should rely on trace 
retrieval algorithms to compute the full autocorrelation trace from a-priori 
knowledge of the electron bunch properties. 



2 Feasibility study 

First we let the electron beam through the first part of the undulator, which 
is 36 m long and is resonant at 0.15 nm. A picture of a single-shot beam 



7 



Table 1 

Parameters for the short pulse mode used in this paper. 

Units Short pulse mode 



Undulator period 


mm 


35.6 


I and II stage length 


m 


35.9 


Segment length 


m 


6.00 


Segments per stage 


- 


6 


K parameter (rms) 


- 


2.9805 




m 


27 


Wavelength 


nm 


0.15 


Energy 


GeV 


17.5 


Charge 


nC 


0.025 


Bunch length (rms) 


/.mi 


1.0 


Normalized emittance 


mm mrad 


0.4 


Energy spread 


MeV 


1.5 



power distribution after the first part of the undulator is showrH] in Fig. 6. 
If one would make an average over many shots, one would obtain (P(t)) 
given in Eq. (1). 

After the first part of the undulator, the electron beam is delayed relatively 
to the photon beam, of a continuously tunable temporal interval t, as illus- 
trated in Fig. 4. Moreover, the microbunching produced in the first part of 
the undulator is washed out. Energy spread and energy loss induced during 
the linear process in the first part of the undulator are taken into account, 
but they are small, and the electron beam can still undergo the SASE pro- 
cess in the second part of the undulator. At this point, the radiation pulse 
P 2 (t, t) is produced. A picture of a single-shot beam power distribution after 
the second part of the undulator is shown in Fig. 7. If one would make an 
average over many shots of this figure, one would obtain (P2(t, x)> given in 
Eq. (2). We performed averaging over 10 shots for each value of x. 

As discussed before, in order to obtain the intensity autocorrelation function, 

3 Note that the Genesis output consist in the total power integrated over full grid up 
to an artificial boundary without any spectral selection. Therefore, Fig. 6 includes 
a relatively large spontaneous emission background, which has a much larger 
spectral width. The second stage of our setup automatically selects the coherent 
contribution within narrow bandwidth only. As a result, the simulation of the 
average power profile at the exit of the first stage can only be performed with the 
help of a special post processing filter. 
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Fig. 6. SASE beam power distribution (single shot) after the first undulator part (36 
m-long). 

after having calculated {Pi{t, t)> for different values of t we need to integrate 
it in time. The result is the average energy in the pulse at the exit of the 
second undulator as a function of t calculated with a 3D FEL code, i.e. the 
autocorrelation trace, which is shown in Fig. 9 (black circles). Averaging 
have been performed over ten shots. 

As shown before, in the ID approximation, the autocorrelation trace cal- 
culated in this way should just coincide with the autocorrelation of the 
average power (i.e. the gain profile) for the first undulator stage. We can 
independently calculate such gain profile. Namely, instead of considering 
the start-up from noise in the first undulator, we can simulate the case when 
a constant laser power is fed at the entrance of the FEL and no shot noise is 
considered. As a result, we obtain (P(t)}, the ensemble average of the power 
distribution after the first undulator for startup from shot noise, i.e. the gain 
envelope. This result is plotted in Fig. 8 (black circles). The expected func- 
tional dependence of (P) on time is given in Eq. (1) and, once normalized to 
unity, it can be equivalently written as 



where I(t) is the current, also normalized to unity for simplicity, I(t) = 
exp[-(t - t ) 2 /(2a 2 )]. Since a and t are known, Eq. (4) can be used to fit the 




(4) 
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Fig. 7. Beam power distribution (single shot) after the second part of the undulator 
(36 m-long)single shot. The SASE seeded power distribution in this plot is obtained 
with a fresh beam and for zero delay. 

simulation data, with a as the only free parameter. It turns out that the best 
fit, shown with a solid line in Fig. 8, is for a = 8.9. 

Since the best fit value for a is now fixed, we can calculate the autocorrelation 
trace using Eq. (3). The result is plotted with a solid line in Fig. 9. It is seen 
that there is a good agreement between actual intensity autocorrelation and 
the black circles. Deviations can be due to differences between the ID and 
the 3D treatments, to the fact that in the second stage we work near to 
the non-linear regime, that we neglect the temporal dependence of Pq, or 
to the fact that for these exemplifications we used only 10 shot averaging. 
However, this accuracy is sufficient for our purpose of demonstrating the 
feasibility of the method. 

Now, let us suppose that we do not know the gain curve, but we simply 
measure the energy per pulse from our setup, i.e. the black circles in Fig. 
9. By inspecting this autocorrelation trace, we conclude that the FWHM 
of the autocorrelation function is about 1.9[im. Assuming, as discussed 
above, a deconvolution factor of 1.5, we obtain an estimate for the FWHM 
of the radiation pulse of about 1.3 fim. The actual FWHM of the average 
power distribution from Fig. 8 is, instead, of 1.4 ^m. This gives an idea 
of the accuracy of the estimation of the radiation profile width with the 
deconvolution factor 1.5. 
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Fig. 8. Normalized gain envelope after the first undulator in the case of startup from 
constant laser power. This plot is identical to the envelope of the average power in 
the case of startup from shot noise. 



3 Simplest measurement using a magnetic chicane only 



As an alternative to the method considered above, we also propose a tech- 
nique to measure the width of ultrafast radiation pulses based on simpler 
hardware. The idea is to rely on a magnetic chicane only, without an optical 
delay line, and is illustrated in Fig. 10. Also a magnetic chicane alone, in 
fact, can provide a delay of the electron beam relative to the radiation pulse, 
as illustrated in Fig. 11. However, the compaction factor of the magnetic 
chicane must always obey to the constraint to be large enough to allow for 
the microbunching produced in the first part of the undulator to be washed 
out. Therefore, the delay t cannot be set to zero as in the previous case, and 
the presence of a simpler hardware is paid by the fact that the setup cannot 
provide a full autocorrelation trace. 

It is therefore necessary to recover the missing data with the help of com- 
puter simulations and information about the electron bunch properties, 
which is available from other measurements. 
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Fig. 9. Energy per pulse recorded at the integrating detector as a function of the 
delay x. It constitutes the intensity autocorrelation trace. 

4 Conclusions 

In this paper we presented a novel method to measure the width of ultra- 
short femtosecond pulses of x-rays produced in XFELs. The technique is 
based on the measurement of the average intensity autocorrelation function 
of the x-ray pulse. To this purpose, we use a fresh-bunch technique. Mea- 
surement of the energy in the radiation pulse produced in our setup yields 
the autocorrelation function. The method works with limited hardware, in 
its simplest form with a weak magnetic chicane installed in place of one 
of the undulator segments. The magnetic chicane can be installed without 
any perturbation of the XFEL focusing structure, and does not disturb the 
baseline undulator mode. Our proposal is therefore cheap, robust, and does 
not present any risk for the functionality of the facility. 
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scheme for X-ray pulse characterization employing magnetic delay only- 
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Fig. 10. Simplest ultra-short X-ray pulse-measurement setup which utilizes magnet- 
ic-chicane delay line only. The setup cannot provide full autocorrelation trace. The 
compaction factor of magnetic chicane must be large enough ( > 1/mi ) for perform- 
ing "fresh" bunch technique. The trace retrieval is based on computer simulations 
and a-priory information about electron bunch properties. 
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Fig. 11. Sketch of X-ray intensity autocorrelation measurement using tunable mag- 
netic delay. This autocorrelator produces part of autocorrelation trace only. 
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